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1 Introduction 



In this paper we consider the aggregation scheme introduced in the paper of Granger 
(1980), where it was shown that aggregation of random parameter AR(1) processes 
can produce long memory. Since this work, a large number of papers were devoted 
to the question how do micro level processes imply the long memory at macro level 
and applications (mainly in economics) (see Haubrich and Lo (2001), Oppenhcim 
and Viano (2004), Zaffaroni (2004), Chong (2006), etc.). 

Recently, Dacunha-Castelle and Oppenhcim (2001) (see also Dacunha-Castelle 
and Fermin (2006)) stated the following problem: which class of long memory 
processes can be obtained by the aggregation of short memory models with random 
coefficients? 

Let us give a precise formulation of the problem when the underlying short 
memory models are described by AR(1) dynamics. Let e = {st, t € Z} be a sequence 
of independent identically distributed (i.i.d.) random variables (r.v.) with Ee t = 0, 
Ee| = <j|, and let a be a random variable supported by (—1, 1) and satisfying 



E 



1-a 2 



<oo. (1.1) 



Consider a sequence of i.i.d. processes Y^ — {Y^,t G Z}, j > 1 defined by the 
random AR(1) dynamics 

Y t u) =aWY<£\+e¥ ) , (1.2) 

where e^) = {e^,t € Z}, j > 1, are independent copies of s = {e t ,t e Z} 

and cS^i's are independent copies of a. Here, the sequences a,a^\ j > 1 and e, 

£ j j > lj are independent. Under these conditions, (|1.2[1 admits a covariance- 
(i) 

stationary solution Y t and the finite dimensional distributions of the process 
= N^ 1 ^ 2 ^2 l I ^ =l Y^ N \ t € Z, weakly converge, as N — > oo, to those of a 
zero mean Gaussian process X t , called the aggregated process (see Oppenheim and 
Viano (2004)). 

Assume that the distribution of r.v. a admits a mixture density ip, which by 
(jl.ip satisfies 

<p(x) 



-i 1 



dx < oo. (1-3) 



The covariance function and spectral density of aggregated process X t are given, 
respectively, by 



j(h) := Cov(X hl X ) = at I j ^(x)dx (1.4) 

-1 J- — 3? 
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and 

Note that an aggregated process X t possess the long memory property (i.e. 

E^Loo \i( h )\ = 00 ) if and onl y if 

f 1 <p(x) 

J_ 1 ^^r dx = °- (L6) 

If the mixture density tp is a priori given and our aim is to characterize the 
properties of the induced aggregated process (moving average representation, spec- 
tral density, covariance function, etc.), we call this problem an aggregation problem. 
And vice versa, if we observe the aggregated process X t with spectral density / 
and we need to find the individual processes (if they exist) of form (|1.2p with some 
mixture density ip, which produce the aggregated process, then we call this problem 
a disaggregation problem. The second problem, which is much harder than the first 
one, is equivalent to the finding of ip such that (|1.5[) (or (|1.4[l ) and (|1.3[) hold. In the 
latter case we say that the mixture density tp is associated with the spectral density 
/■ 

Sections [2] and [3] are devoted to the disaggregation problem. Equality ()1.4|) 
shows that the covariance function "/(h) can be interpreted as a /i-moment of the 
density function ip(x)(l — a; 2 ) -1 supported by (—1, 1), and thus finding of the mix- 
ture density is related to the moments' problem (see Feller (1976)). In Section [21 
we prove that, under some mild conditions, the mixture density associated with the 
product spectral density can be obtained from the "elementary" mixture densities 
associated with the multipliers. In Section [3J we apply the obtained result to the 
spectral density which is a product of to FARIMA-type spectral densities. Using 
the form of mixture density for the classical FARIMA model we show that one 
can solve the disaggregation problem for more complex long memory stationary 
sequences. 

For the aggregation problem, it is important to characterize, for a given class 
of mixture densities, the behavior of the coefficients in its linear representation and 
the behavior of the spectral density. We address this problem to Section 01 In the 
Appendix, we provide a proof of the form of the mixture density in FARIMA case. 
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2 Mixture density for the product of aggregated 
spectral densities 

Let X\ yt and Xi,t are two aggregated processes obtained from the independent 
copies of AR(1) sequences Y\,t — o>{¥\,t-\ + £i,t an d Vz,t = 02^2,4-1 + £2,*, respec- 
tively, where a\, 0,2 satisfy (II. ip . £, = {£j,t, i £ Z} and are independent, i = 1,2. 
Denote of e := Eef t , i = 1,2. 

Assume that c^?i and y>2 are the mixture densities associated with spectral den- 
sities fx and /2, respectively, i.e. 

The following proposition shows that, if y?x an d ^2 in (|2.1[) are supported by [0, 1] 
and [—1,0] respectively, then the stationary Gaussian process with spectral density 
/(A) = A(A)/ 2 (A) can also be obtained by aggregation of the i.i.d. AR(1) processes 
with some mixture density ip and noise sequence e. 

Proposition 2.1 Let tpi and ip2 be the mixture densities associated with spectral 
densities fi and f2, respectively. Assume that supp(<^i) C [0, 1], supp(</?2) C [—1,0] 
and 

/(A) = /i(A)/ 2 (A). (2.2) 
Then the mixture density tp(x), x G [—1, 1] associated with f is given by equality 

v{ x) = ^ (nto ]_ x (1 _ xy){l _ y/x) % + J q {1 _ xy){1 _ y/x) *y 

(2-3) 

where C* := Jq ( f® t tpx(x)(p2(y)(l — xy)~ 1 d]/jdx. The variance of the noise is 

2 2 

CTe = 2n ■ (2 ' 4) 

Proof. Obviously, the covariance function ^y(h) — Cov(X/ l , Xq) has a form 

1 00 

< h ) = V £ 7i(i + N)72(i), (2.5) 

Z7T * — ' 

j=-oo 

where 71 and 72 are the covariance functions of the aggregated processes obtained 
from the mixture densities tpi and ip?, respectively: 

f 1 x^ f° x^ 

71 (j) = °l,e I Pi 0*0 1 2 dx ' ^2 (j) = <J 2 .e / ^(z)- J dx. 
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Clearly, 71 (j) > and 72 (j) = : |7 2 (j)|- Let h>0. Then 

OO 



/ 00 00 ft-— 1 \ 

= IJ l,e cr 2,E ( X! 71 + ^ 72 (•?') + X] 71 ^ 72 + ^) + X 71 (~J + k ^ 2 ^) ) 

S'=o j=o j=i ' 

=: S1+S2 + S3. (2.6) 



We have Si = limN^oo s[ N \ where 



,1,1 1 1 w 

/ / j- ^(-ljV+V'dxdy 

Jo Jo L — x 1 — y 



. , , , x ft 1 l-f-ary)^ 1 , , 



Note that 



1 l-(-*y)"+i 

piOw 2 (-s/> 



1 — a; 2 1 — y 2 I + xy 



< 2<p 1 (x)<p 2 {-y) 



1 — x 2 1 — y 2 



and 



nyi(x)y 2 (-y) 1 1 2 -, 1 2 da; Ay = 71(0)72(0) < 00. 
1 — x 1 — y 

Therefore, by the dominated convergence theorem, we obtain 



f 1 f 1 , * , n X>1 1 i- l-{-xy) N+ ' , , 
S1= L Jo ^ {x) ^ { - V) l-x 2 l-y 2 ^ 1 + xy ^ 

= ll ^ {x)v2{ - V \^ 2 l-y 2 ll X y Axdv 

= J! /_° ^ i-Vi-U d 4 dx - (2j) 



Analogously, we have 



S2 
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For the last term in the decomposition (|2.6p we have 
^ f 1 x- j+h f° xi 

^(fe M-y)j=yz(-y/x) 1 + y/x dx dy 

l fa) f Q rt* h-*i- y /x dx ] dy - (2 - 9) 

Equalities (f2?fj> — ([279 ]) . together with (|2~6| . imply 



/ ^(y) ., 1 2 ^ — dylda; 
J_! l-y'l-xy ) 

y h ( f 1 11 

<P2{y)< / ¥>i(af)i Ji 



-i i - r Uo 1 - x 1 - ^ 

^ X \-x2l-y/x dx } dy 



X n , . T , , 1 







viWi / ^2(^)7; ttt; T d yr da; 



1 — x 2 [ J_ 1 (1 — y/x)(l — xy) 

^2{y){ / ¥>i(af)7i 7T71 rdx}dy. 



l-i I -V 2 [Jo ' - x/y){l- xy) 

This and (11. 4[) imply (|2.3[) . taking into account that J", <^(a;)da; = 1. □ 

3 Seasonal long memory case 

In this section, we apply the obtained result to the spectral densities f± and ji 
having the forms: 

A(A;di) = ^|i-e iA r 2dl (3.1) 

/ 2 (A;rf 2 ) = i |l + e iA r 2d2 , 0<di,da<l/2. (3.2) 

We call these spectral densities (and corresponding processes) fractionally inte- 
grated, FI(g?i), and seasonal fractionally integrated, SFI(e?2), spectral densities. The 
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mixture density associated with the FI(di) spectral density (|3.1[) is given by the 
following expression (for the sketch of proof see Dacunha-Castelle and Oppenhcim 
(2001)): 

¥>i(as;di) = C{d 1 )x^- 1 {\ - z) x - 2dl (l + x)l [0A] (x) (3.3) 

with 

r( » _ r ( 3 ~ d) _ o2d-2 sin(Trd) r(3 - d) 

1 j 2r(d)r(2-2d) r((3/2)-d) 1 j 

and the variance of the noise 

2 sin(7re?i) 



'l,e 



£7(di)7r 



(3.5) 



For convenience, we provide the rigorous proof of this result in Proposition 15.11 of 
Appendix. 

Similarly, the mixture density associated with the spectral density (13. 2|) is given 

by 

^ 2 (a;;d2) = ^i(-a:;d 2 ) = C(d2)|a ; | d2 - 1 (l + x) 1 - 2d2 (l-a ; )l [ _ 1)0] (a ; ) (3.6) 

and 

2 sin(7rd 2 ) . . 

a2 < £ = ccd^T' (37) 

since 



/ 2 (A;d 2 ) = A;d 2 ) = — ^ ]TT ^ ¥ dx = — ^ ^ _ ^ d x . 



Clearly, mixture densities ipi(x;di) and <p 2 (x; d 2 ), given in (|3.3|) . (|3.6p . satisfy 
(jl .3|) and, hence, assumptions of Proposition 12.11 whenever < d% < 1/2 and 
< d 2 < 1/2. Moreover, since di > and d 2 > 0, both <p\ and ip 2 satisfy (|1.6[) . 
The mixture density associated with the spectral density 

/(A;d 1 ,d 2 )=/(A;d 1 )/(A;d 2 ) = -^ |1 - e iA |" 2dl |l + e iX \~ 2d \ 0<d 1 ,d 2 <l/2 

(3.8) 

can be derived from representation (|2.3[) . 



Denote F(a, b, c; x) a hypergeometric function 



F(a,b,c;x) = )' / ^(l - *) c_6_1 (l " 

r(6)r(c- 6) 7 

where c>6>0ifa;<l and, in addition, c— a — b>0ilx = l. 
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Proposition 3.1 The mixture density associated with f(-;di,d 2 ) (|3 - 8[) is given by 
equality 

V {x;d 1 ,d 2 ) = C(d 1 ,d 2 )x d ^ 1 (l~x) 1 - 2dl G(-x;d 2 )l [QA] (x) 

+C{d 2 ,d 1 )\x\ d --\l + xf-^Gix-d^^x), (3.9) 

where 

G(x; d) := F(l, d, 2 - d; ~) - xF(l, d, 2 - d; x) 

and 

TTie variance of the noise is 



o- e = (2tt) G C(di)C(d2)«r 1)E o- 2)E = r-g . (3.10) 



PROOF. (|2.3p implies that 

<p(x\ di,d 2 ) = -J- {C{d 2 )y x {x)F(-x; d 2 ) + C(di)ip 2 (x)F(x; di)), (3.11) 

where 

F(x;d) := r^hOl^d, 

7o -y/x) 

C» = / pifaOf /° dj/)dx = C {d^C (d 2 )C* . 

Jo \J-ii-xy J 

Using equality 

1 + y 1 x 



(1 -xy)(l - y/x) (l-x)(l-y/x) (l-x)(l-xy)' 

we have 

= i f 1 y d - 1 (i-y) 1 - 2d d * f 1 y d - 1 (i-y) l - 2d dy 
l-xj I -y/x l-xj 1-xy 

lunii2 ~ 2 ' ,] 1 (F(l,d,2-d;l/x)-xF(l,d,2-d;x)) 



r(2 - d) 1 - x 
r(d)r(2-2d) G(x;d) 



(3.12) 



r(2 - d) 1 - x 

Now, (|3J)]) follows from (|3TTj) and (f3TT2|) . whereas (|3~T0| follows from (|2^|) . (|3~5 

EH). 
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To finish the proof note that all the hypergeometric functions appearing in the 
form of the mixture density are correctly defined. □ 

In the next proposition we present the asymptotics of ip(x; d\, d 2 ) in the neigh- 
borhoods of and ±1. 

Proposition 3.2 Let the mixture density if be given in (|3.9p . Then 



V {x;d u d*) ~ { C * si « | fd+d 1 (3.13) 

< c ;r 2 &) 1 ~ x ' " (3.i4) 

C* sin(7rdi) 

Proof. Applying identities (see Abramowitz and Stegun (1972)) 

1 T(c)T(c - a-b) 



F(a, b- c; 1/x) = (-^) c - a; c; ), F(o, 6; c; 1) 

\3f — 1/ 1 — X 

we have that for x — > 0+ 

G(-x,d 2 ) = F(l,d2,2-d2; — -) + a:F(l, d 2 , 2 - d 2 ; -2:) 



r(c-a)r(c-6) : 



X 



F(d 2 , l-d 2 ,2- d 2 ; 1/(1 + x)) + xF(l, d 2 , 2 - d 2 \ -x) 



1 + x 

x d2 F(d 2 ,l-d 2 ,2-d 2 ;l) 

r(2-d 2 )r(i-d 2 ) d2 
r(2 - 2d 2 ) x 

^T(2-d 2 ) d2 



2 i-2*r((3/2) - da) 
and similarly for a; — > 0— 



C(r dA ~ V^T(2-di) dl 



This and equality (|3.9[) imply 



y^r(2-d 2 ) _ lli+<Ja _ 1 



V (x;d 1 ,d 2 ) ~ Cfo.fr) 21 _ 2<far((3/2) _ tfa) * 

= g^g^) r(d 2 )r(i-d 2 )x^+^- 1 



sin(7T(i 2 ) 

7T 

sin(7rdx) 



ip{x;d u d 2 ) ~ (C*)- 1 . ^ , ; M^ 2 ' 1 , z^O 
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For x — > 1— we obtain 

<p(x;d u d 2 ) ~ 2C(d 1 ,d 2 ) J F 1 (l,d 2 ,2~d 2 ;~l)(l-x) 1 - 2dl 

ol-2d 2 

sin(7ra2j 

and similarly for x — ► — 1+ 

^(s;di,da) - 2C(d 2 ,d 1 ) J F(l,d 1 ,2-di;-l)(l + x) 1 - 2d2 



Remark 3.1 Clearly, 

^2(2^2 



(^4-^(1 -a:) 1 " 2 * 

sm(7rdi) 



C(di)a; dl - 1 , x->0+, 

2C(d 1 )(l-:r) 1 ~ 2dl , ar^l-, 

C{d 2 )\x\ d *-\ x^O-, 

2C(d 2 )(l + x) 1 - 2d2 , x^-1 



Hence, by Proposition 13. 2i the mixture density <p associated with the product spec- 
tral density (|3.8[) behaves as cpi when x approaches 1, and behaves as tp 2 when x 
approaches — 1. However, at zero, <p behaves as 2_1 , i.e. both densities ipi 

and ip 2 count. 

Proposition 12.11 allows us construct the mixture density also in the case when 
the spectral density / of aggregated process has the form 

1 / IA|\- 2d 
fW = ^ ( 2sin V) 5(A), 0<d<l/2, 

where g(X) is analytic spectral density on [— 7r,7r]. In general, the existence of the 
mixture density associated with any analytic spectral density is not clear. For 
example, AR(1) is aggregated process only if the mixture density is the Dirac delta 
function, what is difficult to apply in practice. Similar inference concerns also the 
ARMA processes, i.e. rational spectral densities. Another class of spectral densities 
obtained by aggregating " non-degenerated" mixture densities is characterized in the 
following proposition. 

Proposition 3.3 A mixture density ip g is associated with some analytic spectral 
density if and only if there exists < a* < 1 such that supp((^) C [—a*, a*]. 
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Proof. For sufficiency, assume that there exists < a* < 1 such that supp(i^) C 
[—a*, a*]. The covariance function of the corresponding process satisfies 

r 1 It|I' 1 I 

\l(h)\ < a\ \ ±±^<p( x )dx 

2 r w w ( , , 



1 - a; 2 

(x) 
1-x 2 



S 2 <P( X ) 

< (T^al / ^ dx 



— a 



= Cal hl , 

i.e. the covariance function decays exponentially to zero. This implies that the 
spectral density /(A) = (27r) _1 J^hL-oo l(h)e is analytic function on [— ir, ir] 
(see, e.g., Bary (1964, p. 80-82)). 

To prove the necessity, assume that / is an analytic function on [— 7r, tt] or, 
equivalently, the corresponding covariance function decays exponentially to zero, 
i.e. there exists 9 € (0, 1) and a constant C > such that \"/(h)\ < C0' h '. Assume 
to the contrary that supp(<ys) = [—1,1]. 

Let h > be an even integer. Then 

h(h)\ = al f z^z<p(x) dx<C6 h 
J —l i — x 

dx<C. (3.15) 



implies that 

\x\\ h <p(x) 



1 — X 2 



Rewrite the last integral as 



V V^) dx = [ 1 \ x \*^Wdx+ f 1 Mfcje^d* 



y/e'-' l-(6x) 2 ~ l-(te) 2 l-(to) 2 



1 - (6x) 



dx =: h{h)+I 2 (h) + I 3 {h). 



For every i£ [ — 1, 1] we have 



1 1 1- (6x) 2 ~ 1- (6>x) 2 ' 

Hence, by the dominated convergence theorem and (| 1 . 3|) . Ii(h) — > as /i — ► oo. 
The Fatou lemma, however, implies that both integrals 12(h) and 13(h) tend to 
infinity as h — ► 00: 

/•V 9 ^(fc) 

liminf^W > / ■: ,„ .„ fiminf x^ dx = 00 

h— oo Ji 1 - (fe) 2 /l-»0O 
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since 1/0 > 1. This contradicts (j3~T5|) . □ 

Example 3.1 Assume that the mixture density <p g has a uniform distribution on 
[a, b], where — 1 < a < b < 1. By Proposition 13. 3[ the associated spectral density is 
analytic function on [— n, n] and can be easily calculated: 

f (x\ a l f b dx 

Jgy ' 2ir(b-a)J a l-2a;cosA + a; 2 

' b — cos An f a — cos A 



fb — cos AN / a — cos ANN 

arctan — — — arctan — — , A j= 0, ±7r 

V sin A / V sin A / / 



2ir(b — a) sin | A| 

/ s (0) = ( r £ 2 (27r)- 1 (l - a)-Hl - b)-\ f g (±n) = ^(2^(1 + a)"^ + 6)" 1 . 

We obtain the following corollary. 

Corollary 3.1 Let ipi(x;d) p.3p and f g (x) be the mixture densities associated 
with spectral densities fi(X;d) (|3.ip and analytic spectral density g(X), respectively. 
Assume that supp(ip g ) C [— a*,0] 7 < a* < 1, and 

1 / |A|N~ 2d 
/(A) = -(2siny) g(X). (3.16) 

Then the mixture density ip(x), x G [—a*, 1] associated with f is given by equality 



where 



dy 



C, := 



1 - ^ / 



4 The structure of aggregated process 

In order to make further inference about the aggregated process X t , e.g., estimation 
of the mixture density, limit theorems, forecasting, etc., it is necessary to investigate 
more precise structure of X t . In particular, it is important to obtain the linear 
(moving average) representation of the aggregated process. 



4.1 Behavior of spectral density of the aggregated process 

In this subsection we will study the behavior of the spectral densities corresponding 
to the general class of semiparametric mixture densities of the form (see Viano and 
Oppenheim (2004), Leipus et al. (2006)) 

ip(x) = (l-x) 1 - 2dl (l + x) 1 - 2d2 il>(x), < dx < 1/2, < d 2 < 1/2, (4.1) 
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where ip{x) is continuous and nonvanishes at the points x — ±1. As it is seen from 
the mixture densities appearing in Section [31 this form is natural, in particular (|4.1[) 
covers the mixture density in (|3.9[) . The corresponding spectral density behaves as 
a long memory spectral density. 

Lemma 4.1 Let the density <p(x) be given in (14. ip . ip(x) is nonnegative function 
on [—1, 1] and continuous at the points x — ±1 with -0(±1) ^ 0. Then the following 
relations for the corresponding spectral density hold: 

/(A) ~ vSLlL) '^ A '" 2 " A -±- ^ 

Proof. Let < A < %. (|4.1|) . (|1.5|) and change of variables u = (x — cos A)/ sin A 
lead to 

fm _ o-l f 1 {l-xf~^{l + xf-^{x) 
/(A) - ^]_ x |l-xe^P dX 



(2 sin-) (2 cos-) 5*(A), 



where 

(T? /-tan A ( gm A u cog A N l-2dW CQS A + ^ gin A N l-2d 2 

5 *(A) = — / 2 2- 2 2- V(wsinA+cos A) du. 

* ./-cot A 1+U 2 

By assumption of continuity at the point 1, the function ip(u sin A+cos A) is bounded 
in some neighbourhood of zero, i.e. ^(usin A + cos A) < Ci(A ) for < A < A and 
Ao > sufficiently small. Hence, 

(sin | -ucos|) 1_2dl (cos| + usin|) 1 - 2d2 u ^ ^ ^ C 2 (A )(1 



V>(w sin A+cos A) < 



1 + u 2 r v ' - 1 

for < A < Ao and, by the dominated convergence theorem, as A — > 0, 

Z7T 

a 2 ^(l) 



2 sin(7T(ii) 

implying (|4.2[) . The same argument leads to relation (|4.3[) . □ 

Remark 4.1 Note, differently from Viano and Oppenheim (2004), we do not re- 
quire the boundedness of function ip{x) on interval [—1, 1]. In fact, tp(x) can have 
singularity points within (—1, 1), see (|3.3[) , (|3.6|l . 
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4.2 Moving average representation of the aggregated process 

Any aggregated process admits an absolutely continuous spectral measure. If, in 
addition, its spectral density, say, /(A) satisfies 



log/(A)dA>-oc, (4.4) 



then the function 



M-)-exp{i-£^±f log/(A)dA}, \z\<l, 

is an outer function from the Hardy space H 2 , does not vanish for \z\ < 1 and 
/(A) = |/i(e lA )| 2 . Then, by the Wold decomposition theorem, corresponding process 
X t is purely nondeterministic and has the MA(oo) representation (see Anderson 
(1971, Ch. 7.6.3)) 

OO 

X t = Y i 4>iZt-i, (4-5) 

3=0 

where the coefficients ipj are defined from the expansion of normalized outer func- 
tion h(z)/h(0), E£o1$ < °°> = !> and Z t = X t - X t , t = 0,1,... (X t is 
the optimal linear predictor of X t ) is the innovation process, which is zero moan, 
uncorrelated, with variance 



a 2 = 2Trexp{i- j log/(A)dA}. 



(4.6) 



By construction, the aggregated processes are Gaussian, implying that the innova- 
tions Zt are i.i.d. N(0, a 2 ) random variables. 

We obtain the following results. 

Proposition 4.1 Let the mixture density ip satisfies the assumptions of 
Lemma \4-l\ Then the aggregated process admits a moving average representation 
(14. 5p . where the Z t are Gaussian i.i.d. random variables with zero mean and vari- 
ance p~eJ|) . 

PROOF. We have to verify that (|4.4|) holds. According to (|4.2|) . log /(A) — 
-Cilog|l - e iA |, |A| 0, where d > 0. For any e > choose < A < tt/3, 
such that 

log /(A) 



Cilog|l-c 



1 > -e, < A < A . 
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Since - log |1 - e iA | > for < A < tt/3, we obtain 

/ °log/(A)dA> (Ci -e) / °log|l-e iA |dA > -oo (4.7) 
Jo Jo 

using the well known fact that J^logll — e lA |dA = 0. Using (|4.3p and the same 
argument, we get 

/•7T 

/ log/(A)dA > -oo. (4.8) 

J it— Ao 

When A E [Ao, 7r — A ], there exist < L\ < L 2 < oo such that 

uniformly in x G (— 1, 1). Thus, by (|1.5p . Li < /(A) < L 2 for any A G [A , n — Ao], 
and therefore 

/ log/(A)dA > -oo. (4.9) 

J\n 



(|47)) - (|4~9]) imply inequality (|4~4)) . □ 

Lemma 4.2 // £/ie spectral density g of the aggregated process X t is analytic func- 
tion, then X t admits representation 

oo 

X t = y^gjZf-j, 

j=o 

where the Z t are Gaussian i.i.d. random variables with zero mean and variance 

a 2 g = 2n^ V [^- J log. 9 (A)dA} (4.10) 



and the 



gj satisfy | J2jL 9j\ < oo, g = 1. 



Proof. From Proposition 13. 31 there exists < a* < 1 such that 



For all x G [—a*, a*] and A G [0, ir] we have 

1 



> C > 0. 



|1 - xe iX \ 2 

This and (|4. 1 1|) imply logg(A)dA > — oo. Finally, IXw^oSjI < 00 f°ll° ws from 
representation 

and the assumption of analyticity of g. □ 
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Proposition 4.2 Let X t be an aggregated process with spectral density 

1 / |AI\- 2<i 
/(A) = — (2sinUj g( X ), 0<d<l/2, 

satisfying the assumptions of Corollaru \3.1\ Then X t admits a representation (14. 5 [I . 
where the Z t are Gaussian i.i.d. random variables with zero mean and variance 

o* = 2*&q>{± y" r iog/(A)dA}=cxp{^ J" log. 9 (A)dA} = gj- 

and the ipj satisfy 



j=o9j .j i , 



(4.12) 



-ffere, t/ie <7j are given in Lemma \4^ 
Proof. We have 



-If 

2vr V 



2 sin ■ 



A| \ - 2d 1 



2tt 



J=0 



with ft. 



r(j + d) 
t(j + i)r(d) 



and, recall, 



E* eliA > E^ 2 <°° 

since, by Lemma |4~21 J* logg(A)dA > — oo. On the other hand, log /(A)dA > 
— oo implies 



3=0 



E^ 2 

J=0 



< oo 



and, by uniqueness of the representation, 



It easy to see that, 



11 7k^' n > !h - 

v 3=0 



E h *-M ~ ftfe E * ~ ^ \ 

J=0 i=0 



(4.13) 



where C2 = V 1 (d)J2'jLo93- Indeed, taking into account that hk ~ T 1 (d)k d l , 
we can write 

& 00 

, 5j 1 : u/;r 'E 

i=o j=o 
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where CJfcj = /ife-^r^fc 1 l{j<fe} — > 1 as k — > oo for each j. On the other hand, 

we have |ofc,j| 5~ C(l + uniformly in k and, since the decay exponentially 
fast, the sum ^°jL Q ()--\- \9j\ converges and the dominated convergence theorem 
applies to obtain (|4. 13|) . 
Hence, we can write 



1P0 



(27T)2| — 

where ipj = 4>jV2~TT / o g ~ C%j d ~ l and 1)4. 5[) follows. □ 



5 Appendix. Mixture density associated with 
FI(ci) spectral density 

Proposition 5.1 Mixture density associated with FI(d) spectral density 

1 / |A|\- 2d 
/(A;d) = — (2si n yj , 0<d<l/2, 

is given by equality 

ip(x) = C(d)x d - l (l - x) l - 2d {\ + x)l [0A] (x), (5.1) 

where 

C(d) = r ( 3 Z d ) = 2 2d-2 sin ( 7rc r ( 3 - d ) 



2T(d)T(2-2d) r ((3/2)-rf)' 
T7ie variance of the noise is 

2 sin(7rd) 
° z = C(d)n ' 

PROOF. Equality 

/(A;d) = 2- d - 1 7r- 1 (l-cosA)-' i 
implies that 1 — cos A = (7r2 d+1 /(A; d))~ 1 / d . Hence, rewriting 

|1 _ xe *\* = (1 - x f(l + (1 - cos A)) , 

and assuming supp(^) = [0, 1], we obtain that the spectral density of aggregated 
process is of the form 



cJS.2) 



2Wo |l-a»iA|a 2^7 (1 - x )a(l + ^^,(1 - cos A)) 

al f 1 <p[x) 



2tt 



o (l-x)a(l+^(7r2*H/(A;d))-V««) 
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The change of variables y x l d — 2x/(l — x) 2 implies 

d2 d x d - 1 (l+x) 

Consider the density tp defined by 



dy = . ^ dXj (5.4) 

C(d)(l-x) 2 V ' 



where C(d) is some constant. Then (|5.2I) becomes 



1 fix) dx = *lC{d) r dy 



2ir J Q |l-a;e iA | 2 2tt J (1 + y 1 /d( 7r 2 d + 1 f(X; d))- l l d ) 



f(X;d)a 2 2 d C(d) 



dz 



l + z i/d 

after the change of variables z = 7r2 d + 1 /(A-d) • Therefore, FI(rf) is an aggregated 



process and, by (|5.3[) - (|5.4[) , the mixture density has a form 

ip(x) = C(d)d2 d x d - 1 {\ - a;) 1 - 2 ^! + x), (5.5) 

and the variance of the noise is (see formula 6.1.17 in Abramowitz and Stegun 
(1972)) 



2 _ _d 



= 2- d (dC(d))- 1 {B(d,l -d)) 
= 2- d (dC(d))-i S -^. 

Finally, it remains to calculate the constant C'(d) to ensure that the mixture 
density ip given in (|5.5[) integrates to one over the interval [0, 1]. We have 

2d dx 



<p(x)dx = C(d)d2 d (j x d - 1 (l-x) 1 ~ 2d dx + J x d {l-x) 1 
= C \d)d2 d i B(d,2 - 2d) + B{d+ 1,2 -2df 
= C(d)d2 1 



^ 1+d r(d)r(2-2d)) 

(2 - d)r(2 - d) 

2 _ d r((3/2)-rf) 



Hence, 



C(d) = 



sm(7ra) 1 (3 — a) 
1 r(3-d) 2 d - 2 sin(7rd) r(3 - d) 



d2 d + l T(d)T{2~2d) d r((3/2) - d) 

and C(d) = C(d)d2 d . □ 
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